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Two- and three-particle distribution amplitudes of heavy pseudoscalar mesons of well-defined geo- 
metric twist are introduced. They are obtained from appropriately parametrized vacuum-to-meson 
matrix elements by applying those twist projectors which determine the enclosed light-cone operators 
of definite geometric twist and, in addition, observing the heavy quark constraint. Comparing these 
distribution amplitudes with the conventional ones of dynamical twist we derive relations between 
them, partially being of Wandzura- Wilczek type; also sum rules of Burkhardt-Cottingham type 
are derived. The derivation is performed for the (double) Mellin moments and then re-summed to 
the non-local distribution amplitudes. Furthermore, a parametrization of vacuum-to-meson matrix 
elements for non-local operators off the light-cone in terms of distribution amplitudes accompanying 
, independent kinematical structures is derived. 
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I. INTRODUCTION 

<D , 

["T i | Exclusive non-leptonic decays of B mesons play a crucial role for our understanding of rare flavour-changing 

processes and the exploration of the mechanism of CP violation within the standard model. Thereby, of 
special interest are hadronic two-body decays, either via B — > D transition, e.g., B — > Dir (DK), or with two 
very energetic light mesons in the final state, e.g., B — > -kit {Kit). While the weak interaction part of these 
processes is fairly well understood their strong interaction dynamics is quite non-trivial. However, some 
simplifications are possible due to the strong ordering of the three fundamental scales, the weak interaction 
0^ ■ scale My/, the 6-quark mass to&, and the QCD scale Aqcd- Because mj 3> Aqcd the heavy quark effective 

theory (HQET) (for a review, see, Refs. 0) may be applied and, furthermore, the strong interaction 
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effects with virtualitics above raj, may be included into the renormalizcd coefficients of local operators Oi of 
the weak effective Hamiltonian. 
\f} , In order to compute the (renormalized) matrix elements {M\M2\Oi\B) , at least in leading order of 

Aqcd/™&, some factorization Q into perturbatively calculable short-distance contributions and appropriate 
long-distance contributions has to be applied - either using the QCZJ-factorization approach |3, |fj |jfl , the 
more effective SCET approach 0,0, or the pQCD approach [Tjj. For example, according to Q that matrix 
element in case of two light mesons can be represented by hard scattering amplitudes T, (BM) form factors 
Ph. Ff M as well as light ($m) and heavy meson light-cone (LC) distribution amplitudes (DA), e.g., 
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(M 1 M 2 \O i \B) =y>f Ml (mI) / duT?,(u)$ Ma (u)+ [ dl; dudvTf'fa u, v) $ B (0 $ Ml («) $m 2 («), (1.1) 
j Jo Jo 

assuming M\ to pick up the spectator quark from the B-meson; obviously, no long-distance interaction takes 
place between M 2 and the (SMi)-system. 

In the case of light (pseudo)scalar and vector mesons the light-cone distribution amplitudes (LCDA) 
are well-known for leading and next-to- lea ding twist for bilocal (quark-antiquark) as well as trilocal (quark- 
gluon-antiquark)operators for the ir- meson |lllll2j and the p-meseon jl3L Il4| - also considering the Wandzura- 
Wilczek (WW) relation. In case of B-mesons they have been determined in the framework of HQET 
also discussing in detail the WW-approximation [l7L[T^ |: the case of the £>-meson is easily obtained, at least 
in leading order, by observing the spin-flavour symmetry of HQET. Furthermore, with the aim of a better 
understanding of the scale dependence of LCDA - and of the hard scattering kernel - in the factorization 
procedure the knowledge of their renormalization behaviour is required. In the case of the leading LCDA 
<1> + this has been studied recently [l^ . 
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In the limit of infinite heavy quark masses (toq — > oo). the heavy quark field Q(x) reduces to an effective 
field h v (x) with the kinematics contained in a phase factor (v being the heavy meson's velocity, v 2 = 1), 

Q(x) -> e- im * vx h v (x). (1.2) 

Moreover, the effective heavy quark field h v (x) obeys the on-shell constraint, 

■fihy = h v , (1.3) 

and the effective Lagrangian £ e g = Y^h^i hvi(vD)hv , is independent of the spin or mass of the heavy quark 
and exhibits, therefore, the U(2Nh) spin-flavor symmetry (Nh is the number of contributing heavy flavors). 
Of course, the on-shell constraint reduces the number of independent heavy meson DAs in comparision with 
light mesons. 

In accordance with the definition of usual meson LCDAs ^lj but additionally respecting the on-shell 
constraint 1)1.3(1 . the B-meson LCDAs arise by parametrizing matrix elements of appropriate non-local LC 
operators Oi which are built up by quark and antiquark fields - occasionally containing also gluons and/or 
quark-antiquark pairs - sandwiched between vacuum (0| and B-meson state \B(v)) of momentum P = M v. 
For instance, the (2-particle) LCDAs are introduced as (phase factor omitted) 

(0\q(K 1 x)rU(Kix,K2x)h v (K 2 x)\B(v))=K r (P,x) J d£ <p B (Q e-^-^H^K (1.4) 

where F denotes some generic Dirac structure 

T = {1,7a: 1^/3,75, 757a: 175^/3} with (J aj3 = ~ [7 , 73] , (1.5) 

and x — x + i>(\/ (vx) 2 — x 2 — (vx)), x 2 = 0, defines some light-ray being related to x by a fixed non- 
null subsidiary four-vector which may be identified with the i?-meson's velocity. The path ordered phase 
factor, U{k\X, K2X) = P exp { ~ig J hl dr x^ 1 A^(tx)} , assuming Schwingcr-Fock gauge, will be omitted in the 
following. Furthermore, the matrix element <|1.4|) is parametrized by a kinematic factor /Cr and the Fourier 
transform of the DA </9_b(£) w.r.t. variable xP; K-r depends on the momentum P of the (pscudoscalar) meson 
and LC coordinate x of the nonlocal operator as well as on the generic Dirac structure T. Explicit forms 
for /Cr will be introduced in Section [H] Everywhere, possible color indices in operator matrix elements will 
be suppressed. The integration range in (|1.4|l results from the fact that, in the framework of non-local LC 
expansion [22j . that matrix element can be shown to be an entire analytic function in the variable xP |23j |. 
Usually, due to the (anti)symmctry of the relevant QCD operators Or{K\X, K2X) w.r.t. exchange ki <-> K2 the 
integration range is restricted to < £ < 1. — Analogous definitions hold in the case of tri- local operators 
including, e.g., the gluon field strength F^^nx) at arbitrary intermediate points kg [k\, K2]. 

Conventionally, LCDAs are characterized by its dynamical twist which, roughly speaking, counts powers 
of M/Q for the various terms in the kinematic decomposition of the matrix elements of non-local QCD 
operators j24|- Alternatively, using group theoretical arguments, the original definition of twist [2^| for local 
QCD operators, r = dimension d— (Lorentz) spin j, has bee n ge neralized to the notion of geometric twist 
for non-local QCD tensor operators on the light-cone [2fil I2H l28| as well as off the light-cone j^E^J. The 
decomposition of such tensor operators into operators of definite geometric twist leads to corresponding 
decompositions of the LCDAs 0, and to their power (or target mass) corrections |3j| • 

In fact, concerning phcnomcnological aspects the notion of dynamical twist is favored. But, from a 
quantum field theoretical point of view, geometrical twist seems to be more appropriate since it has an 
obvious group theoretical meaning and, therefore, should have well-defined renormalization properties; it 
also offers a clear separation between radiative corrections and higher (geometric) twist effects. 

Both definitions of twist, despite of being different for higher twist, coincide at leading twist. However, 
by comparing equivalent kincmatical structures, it has been shown for distribution functions in DIS and of 
LCDAs for light mesons, especially for 7T- and p-mcsons, that there exist unique relationships between the 
distributions of geometrical twist and the usual ones of dynamical twist. 

To be more specific, let us generically denote the distributions of definite geometric twist r by (p\ (rj) 

and the ones of dynamical twist t by Then the distributions of given dynamical twist t are uniquely 

determined by that set of distributions of geometric twist r with r m i n < r < t and, vice versa: 

t T 

^%)= E [drjK/^tpPiv), ^ T \v)=J2 /^(tf-VfoO^tt), (1.6) 

T = T mi „ ^ t=*min 
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with some invcrtiblc kernel Kj l (^,ri), where t m i n = T m j n . In fact, the relations between both kinds of 
distributions are of triangular shape and, therefore, the corresponding set of equations can be solved with 
respect to either basis. Solving the distributions of dynamical twist w.r.t. those of geometric twist allows to 
derive the well-known WW-relations together with additional WW-likc relations |3l| : corresponding relations 
have been derived for the (light) meson LCD As |34| and, later on, called geometrical WW-relations |35|. Using 
in addition the equations of motion, a different set of WW-like relations for the LCDAs - especially for the 
vector meson case - appeared Il4j: thev should be called dynamical WW-relations. This situation has been 
discussed more detailed in Refs. [35j . 

In the case of heavy meson LCDAs the situation suffers from the on-shcll constraint of HQET. First, the 
number of independent LCDAs is reduced. They have been determined already by Grozin and Neubcrt 
[T^ and discussed further by various authors flfil UtI HsI ] , also considering the WW-relation within the 
dynamical twist approach. Second, concerning renormalization, they show some pecularities which have 
been studied more detailed in the case of leading 5-meson LCDA in Refs. [l^- Also regarding this it seems 
to be of interest to consider the geometric twist approach, too, and to find how these two approaches are 
related. The results of such a comparison, mainly based on Ref. |3*3 |. will be presented here. 

The paper is organized as follows. In Sect. ^ we briefly repeat the derivation of the non-local two- and 
three-particle LCDAs of dynamical twist by using the trace formalism according to Refs. 0, fill [2(]]; in 
addition, in order to be able to compare with the corresponding LCDAs of geometric twist, their local form 
are given as Mellin moments. A consistent paramctrization of the relevant matrix elements which is not 
restricted to the light-cone is given in Appendix ^1 hi Sect. IIIII we determine the local two- and three- 
particle LCDAs of geometrical twist by applying the local projection operators (restricted to its light-cone 
form) onto the matrix elements of the corresponding non-local LC operators. In order not to be confused by 
the on-shell constraint and obstacles of renormalization, they are derived for the general case also applying 
to light (pseudo)scalar mesons. The corresponding projection operators are well known from earlier work up 
to tensor operators of second rank [3(1 El; their local form is given in Appendix iBl In Sect. IIVI the relations 
between the LCDAs of definite dynamical and geometrical twist are given, first, for their local form and then 
for their non-local form. Thereby, we also derive the relations between LCDAs of definite geo metric twist 
resulting from the on-shell constraint. Most of the calculations are performed using FORM |37| . 



II. DISTRIBUTION AMPLITUDES OF DYNAMICAL TWIST 



To begin with we briefly review the conventional representation of vacuum-to-meson matrix elements of 
bi- and trilocal light- cone operators for B-mesons in the heavy quark limit by LCDAs of dynamical twist 
times corresponding kinematic structures. Hereby we follow the covariant trace-formalism Q. These LCDAs 
are Fourier transformed and converted into their Mellin moments. 

In the trace-formalism the vacuum-to-meson transition by a generic quark-antiquark operator with a single 
heavy quark is parametrized in terms of two-particle LCDAs &±(vx) as follows |6tll7j: 



(0\q(x)Th v (0)\B(v)) = 



if B M 



iM 75 r 



$ + (vx) 



2(vx) 



(2.1) 



where, as usual, the B-meson decay constant is defined by (0|q , (0)7 Q 75/i„(0)|i3(v)} = i/sMu Q , and M is 
the mass of the B-meson. Computing these traces for the various Dirac structures T one obtains: 



\q(x)j 5 h v (0)\B(v)) = 
\q(xh 5la h v (0)\B(v)) 



if B M | 
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i/flM v a <5> 



\q(x)j 5 ia a ph v (0)\B(v)) = if B M 



(vx) 



2(vx) 



4>_ 



(2.2) 
(2.3) 
(2.4) 



with arguments vx omitted. The parametrization of the matrix element (0\q(x)ia a ph v (0)\B(v)) is 
obtained simply by observing the relation io~ a p = (— i/2)e aj g K yy5i<7 K '\ whereas the matrix elements 
(0\q(x)ia a ph v (0)\B(v)) and (0\q(x)h v (0)\B(v)) vanish. Here, and in the following, we use the notation 



a[ a b/3] ■= (a a b/3 - apba) /2, a^p) ■= {a a bp + apb a ) /2. 



(2.5) 
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As Grozin and Neubert state, in the limit of fast-moving mesons, is of leading (dynamical) twist 
while $_ is sub-leading. When the matrix element l|2.3[) is considered it becomes obvious that v a is related 
to leading twist while x a /(vx) should correspond to 1/Q and therefore is related to subleading twist. 

Equivalently, the vacuum-to-meson matrix element containing a trilocal quark-antiquark-gluon operator 
is parametrized in terms of four three-particle LCDAs ^>a(vx; &v(vx',i!}), Xa^x;^) and Ya(vx;$) with 
■& being restricted to < ■& < 1 as follows [HIU: 



<0|g(i)F^(??i)i ; T/ l , ) (0)|£(?;)) 



Tri 75 r 



l + j) 



(v$ - (ra)7 M ) 



(vx; $) 



- ia^x^^vivx] #) - XuX A {vx;-d) + — — Y A (vx]-d) 

(vx) 



Again, computing these traces for the Dirac structures T (arguments vx and •& omitted) one obtains: 



{{)\q(x)F^(§3:)^ l5 h v (Q)\B(v)) = -if B Mx, 



Xa-Ya 



(0\q(x)F llv (-dx)x' J -f 5 ^ a h v (Q)\B(v)) = \f B M\ (v^Xa - (vx)g m ) ^a - v Xa + 



(vx) 



y a 



(O\q{x)F^£)^ l5 ia a0 h v (O)\B(v)} 



= 2 if B M < g^aXp] ^ v 



(vx) 



^rx [a v fi] Y A + 



(0\q(x)F^($x)x v la h v (0)\B(v)) = - f B M e^x v K x x *y 



(2.6) 

(2.7) 
(2.8) 

(2.9) 
(2.10) 



Again, the parametrization of the matrix element (0\q(x)F^($x)x v \o- a f3h v (Q)\B(v)) follows trivially from 
Eq. (|2.9|l . and {0\q(x)F llv (-dx)x 1 ' 'h v (0)\B(v)} vanishes. According to the above introduced conventions, *$>a, 
and Xa are of leading (dynamical) twist, whereas Ya is of subleading twist. 

An independent derivation of the just introduced parametizations, but not restricted to the light-cone, 
is presented in Appendix ^ In this more general case the matrix elements depend on three leading DAs, 
Xi — X$,X2 and X4, and three subleading ones, Y2, Y4 and Y5, (cf. Eq. (|A.19|) V There it is also shown, that 
off the light-cone "J a and tyy contain also subleading contributions, and <f> contains subleading contributions 
already on the light-conel 

Later on, to obtain relations between DAs of dynamical and geometric twist, the DAs have to be Fourier- 
transformed and, thereafter, be converted into Mcllin moments. As has been argued in Refs. [TH the 
singularities of the DAs in the complex (vx) plane are such that their Fourier transforms &±(u) vanish for 
u < 0. On the other hand, considering non-forward matrix elements of light-cone operators, it has been 
shown with the help of the a-paramctcr representation of Fcynman diagrams, that these matrix elements are 
entire analytic functions with respect to (Px) and, thus, the support of &±(u) is restricted to — 1 < u < +1 

For LCDAs which are related to bilocal operators their Fourier transforms and the corresponding Mellin 
moments read 



$±(wi) 



da 



-iuP. 



(-iPx) r 



<I> 



± In 1 



±\n 



/ du u n $±(u). 
Jo 



(2.11) 



(2.12) 



In case of trilocal operators two-parameter LCDAs occur whose Fourier transforms and the corresponding 
double Mellin moments read 



n=0 



F(vx,#vx) = / Vu t e -i(«i+*«2)^ F( Ui ) = V 
Jo 

n 



(-iPx) r 



m=0 



am p p 
u ± n,mi ± n,r 



T>Ui it" 



m v% F(ui 



(2.13) 
(2.14) 
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here J^Dut = f dui J Q du2 and F(vx; $) generically denotes any three-particle DA. 
After transformation, Eqs. I|2.2jl - (|2.4j) read 

(0\q(x) l5 h v (0)\B(v)} = - ifsMj^ ( ~ lP , 3:)n \ [$+|n + $-|»] , (2-15) 



n 

n=0 



(0\q(x) l5la h v (0)\B(v)) = - if B M ( 1 , " (v a $+\n - rfcr [$+| ro ~ $-|n] ) , (2-16) 



n=0 



n! V" 2(w) 

(-iPS)" U^Zfl 



(0|?(x)7fei* a/J M0)|B(t;)) = i^E" r^T^f [*+|« - *-|»] , (2-17) 

n=0 y ' 

while, the transformed equations, in case of a matrix element of a trilocal operator, read 

(Q\q(x)F^x)x u l5 h v (0)\B(v)) = - if B Mx» [X A |nO?) - Y A \ n (#)] , (2.18) 

n=0 

{0\q{x)F^x)x^ 5 -/ a h v (0)\B(v)) 

= if B Mj2 ( ~ lPa;) " f (« M 5 Q - (ra)^„)* A |„W -5 M ^^A|nW + y^^AlnW) , (2-19) 

n=0 *■ v ' 

(0|g(5i)F AlI ,(^)£ ,y 75i ( T Q/ j^(0)|B(w)) 
°° (— iPa;)" f 

= 2i/ S M^J ^ ^-i {v^V [a X p] + {vx) grfaVfl) [*A|nW - *V|nW] 



n! 

n=0 



n=0 



The local decompositions Ij2.15|l - l|2.21(l w.r.t. dynamical twist are required for comparison with those of 
geometric twist in Sect. II Vl 

III. DISTRIBUTION AMPLITUDES OF GEOMETRIC TWIST 

In this Section we introduce the P-meson LCDAs of definite geometric twist. This is done in analogy to 
the introduction of quark distribution functions in deep inelastic scattering [2?J and of p-meson LCDAs [34| 
in terms of definite geometric twist. Namely, we use the decomposition of non-local light-cone operators 
Os a \ with given tensor structure {a} into a (finite) sum of non-local tensor operators of definite twist r, 

- E °t\ with °t\ = 

T 

with appropriate projection operators, Pj^ CT ^ = vQS a ^(x, d), already known from Refs. [26t l30l l3ll] . 

-gM {<>•'} _ fv(r) ^(r)o' -g(r)[a'/3'] ^(r)(a'/3') \ ,„ „s 

I {a} — y I ' a ' [a/3] ' ' (a/3) ' ' ' ' J ' 

d is the inner derivative on the light-cone (see, Appendix |B|) . Obviously, for a given tensor structure, the sum 
over these projection operators of different twist r defines a decomposition of unity, YliT^ta^" ^ = ^{o-}' 

Considering bilocal operators, the corresponding meson LCDAs of definite geometric twist r, generically 
denoted by <pi T ^(u), are introduced according to Ref. (cf. also Refs. [27l l34| 1 



|O w (i,0)|P( U ))=^pg; {ff ' } (^d)/cW? } («,5)^d u e- iuPi ^\u). 



(3.3) 
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Thereby ICf^f-i (v, x) is the basic kinematical structure (of scale dimension s w.r.t. xd) of that matrix clement 
which can be read off from its parametrization w.r.t. LCDAs of dynamical twist since at leading order 
geometric and dynamical twist coincide by construction. When, in accordance with (|2.12|) . one goes over to 

Mellin moments <f>^\ n one nas ^° a PP^y the corresponding local projection operators V^Jy^^ s (x, d), 

{0 \ O{a/} ( X MB(v)) =Y.f:^:\l:Z& (3-4) 

r n— 

in case of trilocal matrix elements (O|0.{y}(iE, dx, 0)\B(v)) the moment <f>^? n has to be replaced by 

The explicit form of the local projection operators is given in the Appendix [BJ In the following they 
will be applied to the bi- and tri-local vacuum-to-meson matrix elements for pseudoscalar mesons. At first, 
without applying the heavy quark on-shcll constraint. Thereby we obtain a decomposition of these matrix 
elements which differs from the decomposition in terms of dynamical twist determined in the preceding 
Section. Afterwards, we will resum the Mellin moments to get the corresponding i?-meson DAs. 



A. Twist decomposition: Local representation in terms of Mellin moments 

First, we consider the bi- local operators and introduce the Mellin moments of the corresponding DAs 
of definite geometric twist. Then, applying the local twist projectors, in particular (|B.9J) - l|B.ll(l in the 
(axial)vector case and (|B. 12|) - l|B.14|) in the skew tensor case, we get the following decomposition of the 
vacuum-to-meson matrix elements: 

(0\q(x) l5 h v (0)\B(v)) 



(0\q{x)j 5la h v (0)\B(v)) 



(0\q(x) l5 ia af3 h v (0)\B{v)) 



The (pseudo)scalar case is trivial since on the light-cone always p( Ti «"») = 1. Therefore only the single LCDA 

(3) ■ (3) 

rp\n occurs - Obviously, in the axial vector and skew tensor case some LCDAs do not appear, namely, <p A \ n in 
Eq. ()3.6|) and t pf^ n as well as i n Eq. (|3.7|) vanish. The summation in Eq. I|3.7(l begins at n = 1 since ^p^ 
vanishes identically (cf. Eq. I)B.13|) ). The matrix elements of scalar and vector operators vanish completely, 
and that with ia a /3 results trivially from expression l|3.7|) . As stated above, the twist decomposition generates 
all linearly independent kinematic coefficients which arc known from the previous Section, Eqs. - QTty . 
as well as Eqs. I|A.36|) - (|A.38(I of AppcndixE] Let us also remark that, apart from different normalization, 
the leading LCDAs from 1)3. 5f) - 1)3. 7|) coincide with those of Ref . , whereby the additional contribution in 
(|3.6() is considered to be nonlcading (w.r.t. dynamical twist). 

Next, we consider the trilocal operators and introduce the double Mellin moments of corresponding three- 
particle LCDAs of definite geometric twist. Concerning the twist projections, we remind that only the 
tensorial structure of the operators is crucial and not if it is a bilocal or trilocal one. The relevant projection 
operators are defined by l|B.9f) - (|B.11J) for the (axial) vector operators, by 1|B.12|) - 1|B.14(I for the skew 



(r) 



{-iPxT (t) 



r n=0 

oo 



too' {-iPx) n m 



r n—0 

oo 



Tl = 



2{vx) n + 1 



(2) (4) 
f A \n ~ <PA\n 



r n—0 

oo 



(r)[ a 'P'] {-iPx) n yg] ( T ) 

n\ (vx) ^ T l" 



{-\Px) n V [a X p] (3) 

^ n\ (vx) ^ T l™ 

n— 1 v ' 



(3.5) 



(3.6) 



(3.7) 



6 



tensor operators and by (|B.15(1 - (|B.19|I for the symmetric tensor operators of second rank. Unfortunately, 
at present no projectors for tensors of third rank, besides totally symmetric ones, are at our disposal. 
With these operators the trilocal vacuum-to-meson matrix elements in case of a pseudoscalar meson read: 



\q(x)F^(#x)SrT6hv(P)\B(v)) 



>(t)m' 
fi\n-\-l 



if B Mx 



r n—0 

^ (-iPx) n , 



n 



,(5) 



M / y 

n=0 



(3.8) 
(3.9) 



|g(i)F^( 1 ?i)i v 7 57Q /i„(0)|S( 1 ;)) 



n=0 



LPS)" 

— j av*V] 



-CO 
- A1\, 



V 



(na)\n+l n \ 

(-iPx) n 



2(vx) n + 1 



,(4) 
All' 



CO 



(1?) 



-(4) 
- Al|- 



T 



(t) 



t (t) 
1 A2|j 



(0) 



(3.10) 



M a A2\t 



(4) 



-(6) 



|g(i) J F^(7?i)i"7 Q /i„(0)|B(«)) 



/^EEWi 



(t)[mV] (-i-P^)" 



~-fj,' a' kXV % 



n=0 



/bM e„ aK A v k x 



n=l 



(-LPx) n , 



\q{x)F^{dx)x v ^i(j a0 h v {Q)\B{v)) 



2 * M EEn;uim + i 

r ?i— 



(T)(A«'[t»')/8'] (-i-P^)™ / - / ~\ 



(t) 
"T|n 



T£i(0). 



(3.11) 

(3.12) 
(3.13) 

(3.14) 



Here, some remarks are in order: 

First, since the field strength in the trilocal operators 1(3. 8() - (|3.14|) is contracted with x v and, 
consequently, the kinematic terms K, a (v,x), according to relations i(2.18() - ((2.21(1 . have scale dimension 1, 
the local LC projection operators are to be taken for n + 1. Let us remind that in Schwingcr-Fock gauge, 
x v A v (x) = 0, the field strength is related to the gauge potential, nx" ' F^^kx) = — (1 + nd/dnjA^Kx). 

Second, although any tensor of second rank can be split into a symmetric and an antisymmetric part, 
it is impossible to yield an input parametrization for the matrix element in (|3.1U|) with only one set of 
LCD As Y"J"V(i9) associated with the symmetric and another set T^2|„(^) associated with the antisymmetric 
coefficients. The reason is, that the input parametrization has to vanish if additionally contracted with x M . 
Therefore, both LCDAs interfere for the g^a-term. 

Third, as mentioned above, concerning expression 1(3.14(1 we do not know the explicit structure of the 

projection operator ^(L„vg]? n • Furthermore, a projection operator Pn^ii™ does not occur since an 
e-structure on the RHS is forbidden for pseudoscalar mesons. 

Finally, looking at the expressions 1)3. 9|) , (|3.11fl and 13.13(1 we observe again that a huge number of LCDAs 
of definite twist vanishes, thereby having in mind, that an additional free index fj, besides T comes into play. 
In the pseudoscalar case only the highest twist part and in the vector case only T^ n occur (also here 



,(4) 

vto 



0); in the axial vector case, despite of being more complicated, only the LCDAs T 



(4) 



.(6) 



and T „ ., 



occur. Since the on-shcll constraint reduces the number of independent DAs by two, we suppose that in the 
skew tensor case only two additional independent DAs may occur. 

Independently, there occurs another possible set of trilocal vacuum-to-meson matrix elements and their cor- 
responding LCDAs, denoted by f2, which are related to the three-particle operators q(i)F M1/ (-!9i)7 M i l T/i t ,(0). 
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According to their tensor structure they are analogously defined as the bilocal ones, Eqs. i|3.5[) - 1|3.7|> . There- 
fore, we note only their form in terms of Mcllin moments as follows: 



\q{x)F^x)^x v l5 h v {<S]\B{v)) = if B M(vx)J2 



(3.15) 



n=0 



\q{x^F^x)^x u ^ la h v {Q)\B{v)) = if B M(vx)J^ 



(-iPx) n 



n=0 



2(vx) n + 1 



(5) 



\q{x)F^x)^x v l5 ia a0 h v {Q)\B{v)) = 2if B M{vx)J2 



(-iPx) n v [a x 0] (4) 



(vx) 



~\ ll T\n\ V ) ' 



(3.16) 

(3.17) 



Obviously, the matrix elements of equations 1)3. 15|) - 1)3.17)1 are related to those of equations 1)3. 8)1 - (|3.14|l 
by identities of Dirac matrices. In particular, these relations read 



- (0\q(x)F al/ (dS : )x v -f 5 h v (0)\B(v)) 

+ g^{0\q(x)F^(tix)x v 75 ia Xa h v (0)\B(v)), (3.18) 
g x »{Q\q{x)F^{dx)Z v l5lx h v {Q)\B{v)), (3.19) 

\q{x)F^(§x)x»W X h v {Q)\B(y)) 

+ ieAtA a/3 (0\q(x)F^x)x^ x h v (0)\B(v)). (3.20) 

This leads to relations for the corresponding LCDAs which allows to express the O's via the T's. From 
Eqs. (|3.19|l and 13.20fl we receive the following two relations, 



{Q\q{x)F^ v {§x)^x v ^ la h v {Q)\B{v)) = 

(0\q(x)F^(dx)'fx'''r 6 h v (0)\B(P)) = 
(O\q(x)F^x)^x^ 5 ia a0 h v (O)\B(v)) = U a g x p -g%g x a 



,(4) 



(4) 



(4) 



,(4) 



,(4) 
■ A2\: 



(3.21) 
(3.22) 



corresponding to expressions Ti and T2 in Eq. i|A.42|) . However, since we cannot manage the tensor of third 
rank in Eq. (|3.18() . two further relations concerning and f^„($) arc missing. 

The matrix elements containing F fiv ('dx)v fl x l/ follow immediately from the expressions 1)3.8)1 - l|3.14(l by 
multiplication with . Regarding this we should remark that, contrary to the expressions just derived, the 
corresponding result has not the same structure as it would follow from the bilocal operator! 



B. Twist decomposition: Nonlocal representation in terms of distribution amplitudes 

From the local results, we yield the corresponding nonlocal representation by going back to integral 
expressions. The fractions in n transform thereby to a second integral according to 



ib, 



n — r + 1 



f\ uu n f^/vy (3.23) 
Jo J u w \uJ 



Moreover, we have to respect that not all summations include the zcroth moment. Rewriting such sums by 
the exponentials minus the missing moment we use the following formula: 



Lu (e~' iuPi - l)v(«) = J due- iuPi (^b(u) -J ^ S 1>(w) \ . (3.24) 

We thereby yield expressions all multiplied by the same exponential which will be essential, later on, for 
comparison with the distribution amplitudes of dynamical twist. 
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The twist-dccomposcd two-particle distribution amplitudes read in the nonlocal representation 



\q(x)j 5 h v (0)\B(v)) = if B M I d Ui pP(u)e- lu±Jx , 



\q(x)75lM0)\B(v)} = if B M I du <j v a tp%\u) 



2{vx) 



(2) (4) 

<Pa-<Pa 



I 1 - 



w 



(2) (4) 

<P a -<P a 



(w) 



\q(xh 5 ia af3 h v (0)\B(v)) = 2\j B M 



du I flip {u) — 



1 dw ( U \ ( 3 ) 



w \w/ 



—\uP3. 



-\uPx 



(3.25) 

(3.26) 
(3.27) 



(vx) Jo 

The corresponding nonlocal twist-decomposed three-particle distribution amplitudes are given by 

{Q\q{x)F» v {$x)x v l5 h v {Q)\B{v)) = ifsMx^ [ ' Vu^f '(u,,)^" 1 ^" 2 ^ , (3.28) 

Jo 

{0\q(x)F t _ l „('dx)x l, ^ a h v (0)\B(v)) = if B M [ Vui \ (v^x a - g^ a (vx))T^l(ui) - x^VaT^Ui) 



2(vi) 



T (4) T (6) T 4) T (6)l , \ 
1 Al 1 Al A2 ' A2 \ a i) 



x dw 



T (4) T (6) T 4) T (6) 
1 Al ~ 1 Al ~ 1 A2 ~r 1 A2 



(w,u 2 ) ) [ e -^ Ul+€u ^ Pi , (3.29) 



K ~ A 



(0|g(i)F M1/ (^)^ 7ct /i„(0)|B(z;)) = - f B M e^ aKX v*x 



x / VuATf{ Ui ) - I —S( — )T^(w,u 2 ))e 



w \ w 



Ul\ (4). 



-i{u\-\- , du2)Px 



(3.30) 



The nonlocal DAs corresponding to the tensor case (|3.17|) are missing here. 

Since, up to now, we did not apply the on-shell constraint the results H3.25|l - l|3.30|l can be compared 
to the well-known DAs for the 7r-meson |l2j . Of course, up to mass factors which we did not consider, the 
lowest twist DAs coincide. However, due to the factors nj (n + 1) for the moments of higher twist, Wandzura- 
Wilczek-like combinations and, due to missing zeroth moments analogous combinations occur. Furthermore, 
as has been mentioned in the Introduction, the dynamical higher twist contributions are related to geometric 
twist of the same as well as lower order. 



IV. RELATIONS BETWEEN DISTRIBUTION AMPLITUDES OF GEOMETRIC AND 

DYNAMICAL TWIST 

This section is devoted to exhibit the relations between the DAs of definite dynamical twist on the one 
hand and of definite geometric twist on the other hand. Thereby we get also relations among the DAs of 
definite geometric twist due to the heavy quark limit. (Of course, these relations could have been obtained 
by applying the on-shell constraint on the expressions (|3.5|l . I|3.6|l and l|3.7|l in the bilocal case, and (|3.9|) . 
(|3.11(l and (|3.13|) in the trilocal case as well.) By construction, both types of DAs coincide at leading order 
but, in general, differ at higher order. 

First, let us present these relations for the two-particle DAs in terms of Mcllin moments by comparing 
expressions l|2.15ll . (|2.1(i|) and l|2.17|l with expressions (|3.6|l and (|3.7|) . respectively: 

\ ($ +|B + *_,„) = - v f w i (*+,„ - $_,„) = $\ n . (4.2) 

Equations l|4.1|l as well as (|4.2() lead to independent relations for $±. Due to the on-shell constraint in the 
heavy quark limit, the matrix elements of an axial vector operator is related to the matrix elements of a 
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pseudoscalar and a tensor operator. Consequently, we find: 

.(2) _ ,„( 3 ) ,„(3) 



*+|n = = <n-^T„» (4-3) 

*-|« = - ^ + « ^,1) = - (4t + • (4-4) 

Due to these relations, it is sufficient to consider only either the axial vector operator or the pseudoscalar 
and skew-tensor operator without losing information about the heavy quark Mellin moments $±i„. Also the 

relations between the Mellin moments y^ ra , an d <Pt\ti ~ as we ^ as their dependence on <f>±\ n - 

may be read off quite simply. Especially, one finds 



= $ + |o-0>_| , (4.5) 

<P% =-$_|„ + i (*+,„- *_| n ), »>1, (4.6) 

/■i 

^f = -*-|o + y d«ln«($+(«) -*_(«)), n = 0. (4.7) 

Eq. l|4.5(l states $+|o = ^-|oi an< 4 consequently, with (|4.1|l and i|4.2[) . we get ¥>pi = Vaio anc ^ ^t|o = ^- 
Expression (|4.7|) is obtained by taking into account Eqs. I|4.5|) and Ij4.6|) and using l'Hospitals rule. 
By making use of Eq. Ij2.12|l . we are able to obtain the corresponding nonlocal expressions as follows: 

S + (u) = -<p%\u) = V P{u)-<pf{ul (4.8) 

*_(«) =-At.)-/^(A)-A))=-(^) + P f(«)). (4.9) 

The last relation looks like a (geometric) Wandzura-Wilczek relation but with missing non-integrated leading- 
twist term or, stated otherwise, the difference — <!>_ is a pure WW-tcrm in the difference tp^ — ip^ , 

$+(«) - §_(«) = - u 2) ( U ) - ^ 4) ( U )) + (^ 2) ( W ) - ^\w)) . (4.io) 

Equivalently, from the vanishing of their zeroth moments we get some Burkhardt-Cottingham-likc sum rules: 

du $_(«) = f du$+(w), (4.11) 
Jo 

duip { 2\u)= / &uip { p\u), / dutp^'(u) = 0. (4.12) 



/o 

The representation of the DAs of geometric twist by those of dynamic twist reads (omitting trivial ones), 



<pf(u) = §(*+(«) + *-(«)), (4.13) 
^ 3) H = i($ + (u)-$_( U )), (4.14) 



^(u) = -$-{u) + ~J dw($+(w)-®-(w)). (4.15) 

Integrating both sides of (|4.15(l over the range < u < 1 and observing l|4.11(l the result i|4.7|) is re-obtained. 

Now, let us present the relations between three-particle DAs of definite dynamical and geometric twist in 
terms of double Mellin moments by comparing expressions l|2.18|l , (|2.19fl and (|2.21|) with expressions (|3.9|) , 
(|3.11|) and (|3.13|) . respectively. In case of the third order tensor structure we have no projection operator as 
required in the expression H3.14(l at our disposal. Therefore we are unable to present the full set of relations 
between DAs of dynamical and of geometric twist. But, according to the definition of the DAs of geometric 
twist we know at least that relation which results from the identity, <Jr| 5^S^\ of the projection operator 

^Ui[a)i3]\n^ wrncn has to be compared with the corresponding expression in Eq. (|2.20(l . 
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Thereby, we finally obtain the following relations: 



_ T (4) 



*a|„W = T<g |n (tf) = _ T W W + 2(^T) 



.(5) 



T (4) 
A2jn 


T (4) 
X Al|n 


T (6) 
A2|n 


r L Al\n 


r (4) 
1 A2|n " 


T (4) 
1 Al|ra 


T (6) , 
L A2\n 


T (6) 
1 Al|n 



(0) 



(4.16) 
(4.17) 
(4.18) 
(4.19) 



Looking at Eq. H2.20[l we observe that at most two relations connecting \I>y| n (#) and Ya|ti($) with some higher 
twist contributions t!^($), t > 5, are missing. Inverting Eqs. 1)4. 16[) - Q4.19}! the LCD As of geometric twist 
are simply expressed in terms of LCD As of dynamical twist. Thereby, only relations l|4.18|) and (|4.19|) are 
nontrivial leading in the same manner as for the inversion of relation l|4.4|l to vanishing zcroth moments: 



n+1 



Y A \ n {$), n>l. 



(4.20) 



The nonlocal expressions are again obtained by re-converting the double moments to integral expressions. 
Thereby, we use (|2.14f) to get the following nonlocal expressions 

(4.21) 
(4.22) 
(4.23) 



*v(«t) 


= 4 4) K), 






$A(Ui) 


= T^K) 


= T^("i) 


+ T^ 5) ( Ul ), 


X A ( Ui ) 


= T^K) 


= Y A ( Ui )- 




Y A { Ui ) 


= T^K) 


f T^ 5) K) = 


1 r T (4) 

2 [ ± A2|n 



,(4) 
■All 



-(6) 
- A2I 



, T (6) 



1 ^dw I" (4 ) _ (4 ) _ ( 6 ) „(6) 

2 / „, A2|n X Al|n L A2\n ~<~ 1 Al\n 
Jui w 

with Y A (ui) being a pure WW-term. Omitting the trivial relations, the inversions read, 

nl du> 



|(T^( Mj )-T^K)) = \{^ A (u i )-X A (u i ))+Y A {u i ) + j ^-Y A (w,u 2 ), 

rf( Ui ) = Y A ( Ui ) - X A ( Ui ), 

r^\ Ul ) = ^ A ( Ui )-^v(ui). 
In the same manner as for the bilocal case we obtain Burkhardt-Cottingham-like sum rules, 



/ Y A { 
Jo 



Ui) = . 



= 0. 



In addition, due to vanishing of Sl$j {&) = 2T^ (i?) + T^O?) + T^ |0 (tf) and T^tf), we get 



/■l 

/ P Uj T^u,) = 0. 
Jo 



Pui fi^ 4) (uj) = 0, 



"^(t^uO + t^K)) 



the last relation is consistent with expressions 1)3. 10|) and (|3.11|) . From these relations we obtain 



Vui^ v {ui) = 0, f Vu t {^ A { Ui )+ X A {u t )^ = 0. 



(4.24) 

(4.25) 

(4.26) 
(4.27) 

(4.28) 

(4.29) 
(4.30) 



Furthermore, in complete analogy to (|4.7|l . from relation 14.20(1 together with 14.28fl and l|4.30|) we get, 

ijf Vui fr^liu^-T^liui)) = J Vu l V A (u i ) + J Vu i lnu 1 Y A (u 1 ,u 2 ). (4.31) 



For the remaining LCD As Q^p\ui) , 0,^' {Ui) and Q^'iui), because of their similarity with the bilocal 
LCD As, one gets relations analogous to Eqs. (|4.1|) - (|4!l5|l with the <p's replaced by f2's (at higher twist). 



)(3) 



(5), 



11 



V. CONCLUSION 



Two- and threc-particlc LCD As of definite geometric twist have been introduced and discussed in the case 
of B-mesons. Comparing them with the corresponding LCD As of dynamical twist |T^ we were able to derive 
relations between these different types of DAs as well as those relations which are due to the heavy quark 
limit halving the number of independent DAs. Some of these new relations, especially those for <t> + — $_ and 
Ya are of pure Wandzura-Wilczek type. In addition, since some zeroth Mcllin moments vanish various sum 
rules of Burkhardt-Cottingham type appeared. Concerning conventional LCDAs, this are the expressions 
(14.1111 . 14.28(1 and ((4.30(1 which, to our knowledge, up to now have not been considered in the literature. 

In principle, applying the heavy quark symmetry analogous relations for heavy vector mesons could be 
obtained. However, these relations would be more complicated than those considered here. Concerning 
the two-particle LCDAs they could be obtained from the already known relations between geometric and 
dynamic twist for the /9-meson [34^ by applying the on-shell constraint . 

The derivation of these new DAs makes use of projection operators onto local light-cone operators of 
definite twist, sandwiched between vacuum-to-meson matrix elements. Therefore, it was necessary to work 
with (double) Mcllin moments and re-sum afterwards into nonlocal LCDAs. The twist-decomposed three- 
particle LCDAs for third rank tensors are missing since the corresponding twist projectors were not available. 

The projection operators onto local operators of definite twist, given here on the light-cone, are also known 
off the light-cone |3J|. This opens the possibility to consider, at least in principle, mass corrections analogous 
to an earlier study for p-meson DAs Of course, such corrections will interfere with and, therefore, 

supplement the mass expansion in HQET. Concerning B-mcsons such corrections could be obtained, at least 
partly, from an earlier consideration of 7r- meson DAs by Lazar ( |28l| . p. 68). 

We expect that the discussion of the relations between LCDAs of dynamic and geometric twist improves 
the understanding of B-mcsons. In addition, since the LCDAs of geometric twist behave well-defined under 
renormalization, this approach may help to extend the study of renormalization properties of leading twist 
LCDA ll<| also to the non-leading and three-particle ones. 

In a subsequent paper, we are going to establish relations between two- and three-particle DAs forced by 
the equation of motion. In order to do so the knowledge of two-particle DAs off the light-cone is required. 
A first step in that direction has been made in Appendix A where the parametrization of vacuum-to-meson 
matrix elements is considered off-cone. 
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APPENDIX A: DERIVATION OF INDEPENDENT DISTRIBUTION AMPLITUDES OF 
DYNAMICAL TWIST IN THE HEAVY QUARK LIMIT 

In this Appendix we derive, in the heavy quark limit, the linearly independent kinematic structures - 
together with Lorentz invariant DAs in x-space - which are compatible with the tensor structure of the 
various bi- and tri-local light-ray operators when sandwiched between vacuum and S-meson state. Thereby, 
for the sake of convenience, all the Dirac structures ((1.5(1 are taken into account, considering ^\<J a fi and its 
dual ia a p = (—i/2)e a p K \-f 5 i<j KX on equal footing. Due to the on-shcll constraint 1(1.3(1 and because of the 
well-known relations on Dirac matrices we are enabled to derive all these parametrizations from a single 
general ansatz. By flavor symmetry, that parametrization is valid for any heavy meson. In addition, since 
in HQET spin symmetry relates pseudoscalar and vector mesons, a single ansatz is sufficient for both. 

Furthermore, the general structure of matrix elements for bilocal operators (0\q(x)Th v (0)\B(v)) may be 
obtained from the trilocal operators (0\q(x)F fJi ^(dx)Th v (0)\B(v)) by multiplying it with v^x v / (vx), ignoring 
the dependence on {fx and renaming the obtained DAs. Of course, also matrix elements built with the dual 
field strength can be obtained this way. The consideration of trilocal operators is necessary for the study 
of Wandzura-Wilczek relations in the framework of dynamical twist. For the same reason it is necessary 
to know as much as possible about the kinematic structure and the DAs off the light-cone. Therefore, we 
restrict our consideration to the light-cone only at the end. 



12 



Let us begin with the trilocal operators. In order to be able to apply the on-shcll constraint ((1.3(1 we choose 
to parametrize the following matrix element: 

(0\q{x)F liv {'dx)^ p \a al 3h v (Q)\B(v)) = (i/ B M) (F^j 5 j p icr a f3). 

Here and in the following, matrix elements are simply represented by the different parts of their operators, 
also omitting the various variables. The most general ansatz for the independent kinematic structures which 
can be built for a tensor of rank five with two pairs of antisymmetric indices reads 



(Fpv 757pi (J a/3} = 9a[p 9v\f) 9po Z° + 9a{p. 9u][a 9(3]p Z 2 + 9a[a 9(3][p 9u]p Z% 

~ux~ 9p[a 9(i]a 4 ~^x~ 5p[M 9v]a 5 (vxf pa 10 

f 1)^1) X^ 1 X^ 'D U^ X^ X^ \ 

+ 9,[p 9 V]T 9x la 9p] T 9 P * \-jr Z% + — Z1 + -— Zi + , (A.l) 

where Zf := v a Xi + (x cr /vx)Y t ,i = 1, ...,10, with X { = Xi(vx, v 2 , x 2 ; t?) and Y { = Y l (vx,v 2 , x 2 ;^) are 
altogether twenty linearly independent, Lorcntz- invariant three-particle amplitudes of equal dimension and 
parity; furthermore, in this Appendix, the bracket notation 1)2.5(1 is used without the factor 1/2. 

Due to the on-shcll constraint ((1.3(1 and the well-known identities of gamma matrices, the various matrix 
elements for different Dirac structures can be related to (|A.lll . In particular, it holds 



(F^jpio-ap) 


= (l/2)£aPaT(Ff J ,vJ 5 J p io- CrT ), 


(A.2) 




= (1/6)6^(^7^0, 


(A.3) 


<^75> 


= v a (F^j 5 j a ), 


(A.4) 


{Fp V ^ b ia al3 ) 


= v p (Fp^-fpiaap) + 2 (v a (F^ v 757a ) - vp{F llv ^ a )) , 


(A.5) 


{Fpvla) 


= (1/6)6^(^757"^), 


(A.6) 


(Fp„) 


= v a (Fp„ la ), 


(A.7) 


{FftvicTap) 


= v p (Fp U ^pio- a p) + 2 (v a (F^jp) - vp(Fp^ a )) . 


(A.8) 



Relation (|A.2J> avoids the introduction of dual amplitudes and by equations (|A.3(I - (|A.8J) the parametrizations 
for all the basic Dirac structures 1(1.5(1 can be derived from the general ansatz ((A.l(l . 

Up to now, these parametrizations do not respect those dependencies due to the identities for the gamma 
matrices. In fact, it suffices to require the Chisholm identity, 

7^7° = {9 pv 9 a!i ~ 9 pa 9 V& + 9 pf> 9 ua ) 7/3 + ie pvafi l5l p , 

for the LHS of fPjl . 

{Fpulhlp^al3) = 9l3p(Fp^5la) - 9ap(Fp^5ll3) - ^paPa(Fp^ a ) . (A. 9) 

This requirement results in a linear system of algebraic equations for the DAs Xi, Yi with unique solution: 
A3 = X±; Y3 = Yi; Yq = —X-; = X5, Y$ — ~Xg = Y§, 

A 6 -r 7 = A 8 = y 9 = A 10 = r 10 = o ; (A.io) 

whereas the remaining eight DAs X\, Y\, X2, Y2, X4, Y4, X5, I5 are independent ones. 

Applying these restrictions to ((A. 1(1 . computing l(A.2() and the six matrix elements l(A.3() - (IA.8I) one gets: 

(Ffiv'Ts) = ~ v [p9v]<jZ2 - — v [ti x v] v o Z l (A- 11 ) 

{Fpulbla) = - g a {ti9 v \<7 Z 2 ~ -^v^x^g aa Zl (A.12) 
{Fp^ 5 \a a p) = v p e pa p T e T fll/a Z° - v^g^g^ Z°~ 

— -V[pX u] v [a g p]a Zl + —v K -x x v p e pal3T e T KX[ pg„ ]a Z£ (A.13) 
ux ux 

(^^757pi f7 a/3) = e pa f3 T e T p l/a Z{ + g p [ a gf3][fj,g v ]a Z. 



2 

+—V[pX u ]g p [ a g l3 ] (J Z% + — v K x x e paf3T e T K x[p9u]a %l ( A - 14 ) 
UX ux 
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H v^x^iepapa Z1 v h, x x g p [ a ie^ K x[p9iy]cr Z\ 



<>.)' 



(Fn V ia a p) = V[ a iefl MVIT Z\ + iv p e pa p\ ll g v ] (T Z% 



+— V[ tl x v \\v p e pa p a Z% + — v K x A iv[ a e ] K x[^g^a Z$ 
ux ux 



(Fttvla) = '^nvourZl + —iga[nev]a K \V K X Zg 

(Fpu) — itpupcjV^ Z-^ , 



(A.15) 

(A.16) 

(A.17) 
(A.18) 



i.e., one gets eight independent distribution amplitudes. Since v is dimensionless, all the matrix elements 
are of (mass) dimension M 3 - two dimensions are due to Fp V and one due to the bound state and, therefore, 
the structure functions are ~ M 2 . Furthermore, concerning dynamical twist, the DAs Yi are subleading in 
comparison with the DAs X{. 

Now, let us truncate expressions (|A.11(I - (|A.18(I with x v . Thereby, Y\ disappears and, effectively, only 
the following four combinations remain as independent ones (tt = v 2 x 2 /(vx) 2 ): 



*y = X 1 - X 5 - QY 5 , V A =X 2 -nY 4 , X A = -X 4 , Y A = Y 2 + Y 4 . 
As result we obtain: 



(Ff,vx"7 5 ) 
{F^x" 757 Q ) 
(Fp, v x u j 5 ia a i3} 



( Xm ~ h v ") ( Ya ~ Xa ^ 

(vpX a - (vx)g fia ) ^ A + — {x^Xa - x 2 gp a ) Y A - (x^ - 
(vp,X[ a - (vx)g lt[a ) vp] (*y - *a) - — (xpx [a - x 2 g ft { a ) v 0] Y A + g^x^ *y (A.22) 



v. v a X A 



(A.19) 

(A.20) 
(A.21) 



(Fp V x^ 5 jpia a0 ) = 5 p[q 5/3] a {x x v f _ t - (vx)g x ) V A + g p[a gp] x (x x x f _ t ~ x 2 g x ) — Y A 



-e pa pre T p K \v K x *y - g p [ a gp] K v K yxp - —v^jX A 
{Fp U x v "f p i(j al 3) = - ie pQ/ 3A {x x v ll - (vx)g x ) ^ A - ie pal 3\ {x x x l _ L - x 2 g x ) ^- Y A 

2 

+ig P { a £f3] f iK\v K x x *y + ie pa p K v K (x^ - ^v^X A 



(A.23) 



(A.24) 



A 2 A 

^)(*v-*a) + XX "~ X9 ' 1 Y 



(F^x"i<r a p) = ie a p K xv K (x X Vp - (vx)g x ) (*y - * A ) + 
(F^Ja) = ie fiaKX v K x x *y 



v.v 



(F^x") 







itp a fi\x x *y (A.25) 

(A.26) 
(A.27) 



Obviously, according to the antisymmetry of F^ v , after a further truncation with x^, any expression vanishes. 
In the course of the computation we have used 

ie al3KX v K (x x Vp - g x {vx)) - ie^ a px x v 2 = - ie a p KX (g*g x g p + g^g x g p + g^g x g p }jx u Vpv' J , 
{v^x [a - [vx)gp [a ] vp] + v 2 g fl[a xp ] = g a[K gx]p (g^g x g p + g»g x gp[ + g^g^gt) x"v p v a , 

together with the well-known relations concerning products of e-tensors and its contractions, especially, 

£pa/3T€ 7 fiua = 9po-fl'a[/ifl , i/] l 9 + Sp^Si/] [a 0/3] it ■ 
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Now, truncating also with w M the matrix elements <|A.20f) - l|A.26() simplify considerably: 



(i^xVui) = - (vx) (1 - n) $ 2 , (A.28) 

(F^V^a) = a:a*i-«a(«a0[*i + (l-ft)*3], (A.29) 

(F^^gia^) = (Vaxp-vpxa)*! = - (i/2) e a/3reX {F^tW*} , (A.30) 

(F^^ 7c( ) = = (F )MV x v v fJ ') , (A.31) 

with $i = + Ya, $2 = -3d ^ Ya; truncating with g MQ , we obtain 

(F^x»r l5 ) = (vx)(r 2 + nr 3 ), (A.32) 

(F^z<V757a) = z« Tj - Ua («a;) [Ti - (T 2 + QT 3 )] , (A.33) 

(i^V75i<Ta/3> = (« a aJ/J - «/ja:a)Ti = -(i/2)e aM (VY^ A ), (A.34) 

(F^VTa) = = (F^V), (A.35) 

with Ti = 2* y + * A + X A , T 2 = 3* A + X A , T 3 = 3Y A - X A . 



Despite of the different definitions of the DAs <6j, % = 1, 2 and T,-, i = 1, 2, 3, we observe a striking similarity 
in the structure of the two kinds of matrix elements (|A.28(1 - (|A.30|) and (|A.32I) (|A.34|) . 

Let us now consider the case of bilocal operators. The matrix elements (r) can be read off from the 
relations l)A.28|) - l|A.31|) . taking advantage of the similar tensor structures and assuming the two-particle 



DAs not to depend on d. After division by vx one gets 

( 75 > = -($+ + $_)/2, (A.36) 

(757a) = -v a <S>+ + ($+-*_) , (A.37) 

(75i^) = Va \ v ^ X<X (*+-*-) = -m)e a0KX (ia^), (A.38) 

( 7q ) = = (1), (A.39) 



where, for the sake of convenience, we introduced $>±(vx, v 2 , x 2 ) = ((1 — il)<I> 2 ± $i)(vx,v 2 ,x 2 ), i.e., the 
additional x 2 -dependence from £1 = v 2 x 2 / (vx) 2 has been included into the definition of the new DAs. 

Now, let us restrict onto the light-cone, x 2 — > x 2 = 0. Then the linearly independent LCDAs, which do not 
depend on x 2 , will be written with a "hat", i.e., X(vx]i)), Y(v£ti; and so on. The kinematical coefficients 
of relations (|A.20() - (|A.27|) are shortened, and the distribution amplitudes are reduced as follows: 



* V =X 1 -X n , $ A = X 2 , X A = -X A , Y A = Y 2 +%, (A.40) 

$ 1 =*a + Ya, $2 = X a -Y a , (A.41) 

f i = 2$ v + if A + X A , T 2 = 3*a + A > a , t 3 = 3Y A -X A , (A.42) 

$+=X 2 -X 4 , = X 2 + X 4 + 2(Y 2 + Y 4 ) . (A.43) 



Making use of these conventions and the restriction to the LC, the above derived re pres entations of the 
various matrix elements completely coincide with those obtained by Kawamura et al. |l7J using the trace 
formalism, cf. Eqs. (|2.7|) - l|2.1L)[l and (|2.2f> - (|2.4|l in the trilocal and bilocal case, respectively (cf. Subsection 
ITT}) . Hence, we named the invariant DAs already in accordance with that Reference. 

APPENDIX B: ON-CONE PROJECTION OPERATORS ONTO GEOMETRIC TWIST 

In Sect. Hill we made use of (local) projections onto LC tensor operators of well-defined geometric twist. 
The general procedure of decomposing non-local QCD tensor operators, either on-cone or off-cone, into a sum 
of such operators having definite geometric twist has been developed in a series of papers |2a . 1271 \2&L l29l l30| . 
There, it has been shown that this twist decomposition crucially depends on the tensorial structure of 
the operator under consideration and that it can be obtained by using appropriate projection operators. 
That procedure makes use of the representation of non-local tensor operators into a series of local ones and 
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the decomposition of local tensor operators w.r.t. irreducible representations of the Lorentz group. The 
procedure simplifies if light cone operators are under consideration. Since local and nonlocal LC expansion 
are related mutually we can freely choose if we calculate the twist decomposition in the local or the nonlocal 
representation. 

The local LC tensors which are to be decomposed into tensors of definite twist are given, according to the 
representations l|H.4[) . by /Cr S J.M (v, x) (Px) n . In the following they are denoted as 

N n (x), O a \ n (x), M[ a p]\ n (x), M( Q/3 )|„(a;). (B.l) 

in the (pseudo)scalar, (axial) vector, antisymmetric and symmetric tensor case, respectively. Below, we 
present the corresponding local LC projection operators. The LC projection operators for tensors of third 
rank are available only in the totally symmetric case, but this is of no relevance here. 

The most compact representation of the LC projections makes use of the 'interior' derivative |3Sl ] acting 
on the light-cone, 

d a = {(l + (xd))d a -±x a u} , (B.2) 

X— X 

which, together with x a ,X := 1 + (xd) and X a p := xpd a — x a dp spans the conformal algebra so(4, 2). In 
that order, these operators are the generators of special conformal transformations, translations, dilations 
and rotations, respectively. Especially, it holds 

[d a ,xp] = g a pX + X a(3 , (B.3) 

[X, X a ] = X a , [XaffjXfi] = g^aXf) - g^pXa, (B.4) 

[X, da] = - d a , [X a p, d M ] = g^adp - g^pda . (B.5) 

From this it follows 

{X-l)d [a x 0] =-{X + l)x [a d 0] , (B.6) 

d( a xp) = x^dp) + Xg a/3 . (B.7) 

The various LCD As of definite geometric twist as well as their moments are labeled by r = tq + r (2^| . 
Thereby, To is that part of the twist which corresponds (or would correspond) to the totally symmetric tensor 
operator, and r = 0, 1, . . . , labels higher order contributions due to the actual symmetry type characterizing 
the irreducible representations of the orthochronous Lorentz group which appear in the decomposition of the 
light-cone operators. In fact, for the operators ipi(x)Tip2(Q) with T = l,^ a ,\a a p we obtain tq = 3,2, 1 + 1, 
respectively, namely, if the minimal twist of M^ a ^ is To then the minimal twist of M[ a p] is to + 1. In principle, 
there may occur different LCD As of the same twist t accompanying equal kinematical structures. This takes 
place for tensors of higher rank but will not be made explicit here, i.e., only their sum will be given. 

Now, we state all generic non-vanishing local on-cone operators of definite geometric twist up to tensors 
of second rank [30, |3l| : 

N^(x) = N n {x) (B.8) 
°i|n +0) ( 5 ) = J^TjS^O^ix), (B.9) 

= (5 - (^TTF + ia *")) (R10) 



o 



i"^) = Jn+Tp ^O^x), (n > 1), (B.ll) 

= - {n+ i) {n+2 ) ~ d ^ { s m jdw + ^)) ^ Im [hi«(^). (B.12) 

+ (7i+l) 2 (n+2) a[Q<5 S rl + rt^*l<* S 0]* U] ) M ^\ n W' ( n ^ X )' ( B ' 13 ) 

<S, 3 '(-) = -^klfl* - h + ^)) & M M|n(S), (n > 1), (B.14) 
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M £Sn(*) = {n+l)kn+2Y d a d x»x»M { , u){n (x), (B.15) 



M 



= ^fcl) d (« (4) - (n+2) 2 ( d 0)* (M + ^)d ( ^)) ^M^iJx), (n > 1), (B.16) 
SSSC*) = + n(n+l) 2 (n+2) (d«d^ + ^dj - ^ d (a <^ 



~ (n+l)(n+2) x («°/3) u + n 2 (n+2) 2 x (« u /3) x u + n(n+2) 2 x (« u /3) c 

+ + V'"' V(/«0I»(*)> (« > 1), (B.17) 

M &)|n( f ) = (n+lKn+2) £ <<" (4) " ( d « i( " + x «d (AI )) ,„(*), (n > 1), (B.18) 

Af&tlS^) = n 2( n 1 + i)2 £a£/9d"d*'M (/ll , ) | ft (x) ) (n > 2). (B.19) 



Let us remark that the restrictions in n, appearing in Eqs. (|B. 13|> and i|B.14|) as well as Eqs. i|B.16|) - (|B.19ll 

are automatically fulfilled due to the definitions of these expressions, i.e., the zeroth and first moments, 
respectively, vanish by construction. Furthermore, the second term in the RHS of Eq. I)B.12|) and the last 
term of Eq. (|B.14|) vanish - they are written only because of analogous terms in case of the symmetric 

tensor. Finally, we should mention that M^pyZ* (x) contains two and M^ ^^(x) contains five independent 
components corresponding to irreducible representations of the Lorentz group which, in principle, could be 
accompanied by independent DAs. But in that paper we associate only one and the same with them. Any 
other expression corresponds only to a single irreducible representation. 

Obviously, these operators of definite twist are obtained by applying the corresponding projection operators 
Vn , 'PaVn' ^[ajittn > ^{a^ln (including the fractions in n) on the undccomposcd operators: 

= (v ( n T) N n ) (x), (B.20) 

( Z(S) = (rit l O, ln ) (x), (B.21) 

<V*( £ ) = (*). (B-22) 

M &)\n® = (r$ffMQ»>)\n) (*)■ (B.23) 
In addition, they obey the common property of projections: 

(pM x V^)™ = 5 TT V ( p^' n ', (B.24) 

T max 

]T V {t) = 1 (B.25) 



In order to prove these properties the conformal algebra and the relations (|B.6(I and IB. 7(1 have to be used 
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